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Abstract 

Using pp collision data collected by LHCb at center-of-mass energies ^/s = 7TeV 
and 8TeV, corresponding to an integrated luminosity of 3fb~^, the ratio of the 
branching fraction of the —)• ?/)(25)7r+ decay relative to that of the —)■ 

decay is measured to be 0.268 ± 0.032 (stat) ± 0.007 (syst) ± 0.006 (BF). The first 
uncertainty is statistical, the second is systematic, and the third is due to the 
uncertainties on the branching fractions of the J/ij: —and '4>{2S) —)• 
decays. This measurement is consistent with the previous LHCb result, and the 
statistical uncertainty is halved. 
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In the Standard Model of particle physics the meson family is unique because it contains 
two different heavy flavor quarks, charm and beauty. The ground state of the Be meson 
family has a rich set of decay modes since either constituent quark can decay with the 
other as a spectator, or they can annihilate to a virtual W boson. The search for new B'^ 
decay channel^ and precise measurements of their branching fractions can improve the 
understanding of Quantum Chromodynamics (QCD) and can test various effective models. 
Many properties of the B'^ meson have been investigated by the LHCb experiment; the 
Be mass, lifetime and production rate have been measured while several new decay 

The observation of the —>■ '0(2S')7r+ decay 
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channels have been observed 
was made with pp collision data at a center-of-mass energy of -y/s = 7 TeV, corresponding 
to an integrated luminosity of 1.0 fb“^ j^. The ratio of the branching fraction of the 
B^ —>■ '0(2S')7r+ decay with respect to that of the B^ —)■ J/'07r+ decay, defined as 


_ g(i?+ ^ ^(2^)7r+) 
B{B+^ J/'iljn+) ' 


( 1 ) 


was measured to be 0.250 ± 0.068 (stat) ± 0.014 (syst) ± 0.006 (BF). The hrst uncertainty 
is statistical, the second is systematic, and the third is due to the uncertainties on the 
branching fractions of the J/V' —?• and 'ip(2S) —)■ p~^pL~ decays. The statistical 

uncertainty is dominant. Several theoretical predictions for based on different effective 


models 14-19 exist, and vary between 0.07 and 0.29. 


The analysis presented here updates the previous LHCb measurement of Rb using 
the full pp collision data collected by LHCb in 2011 and 2012 at ^/s = 7 TeV and 8 TeV 
respectively, corresponding to an integrated luminosity of 3fb~^. Due to the increased 
data sample and an improved analysis method, the statistical uncertainty is reduced by 
half, allowing a more powerful test of the theories. 

The LHCb detector 20, 21 is a single-arm forward spectrometer covering the 

5, and is designed for the study of particles contain¬ 


LHCb detector 
pseudorapidity range 2 < 


p < 


ing 6 or c quarks. The detector includes a high-precision tracking system consisting of 
a silicon-strip vertex detector surrounding the pp interaction region [^, a large-area 
silicon-strip detector located upstream of a dipole magnet with a bending power of about 
4Tm, and three stations of silicon-strip detectors and straw drift tubes 23 placed down¬ 
stream of the magnet. The tracking system provides a measurement of momentum, p, 
of charged particles with a relative uncertainty that varies from 0.5% at low momentum 
to 1.0% at 200GeV/c. The minimum distance of a track to a primary vertex (PV), the 
impact parameter, is measured with a resolution of (15 -1- 29/pt) where px is the 
component of the track momentum transverse to the beam, in GeV/c. Different types of 
charged hadrons are distinguished using information from two ring-imaging Cherenkov 
detectors |^. Photons, electrons and hadrons are identified by a calorimeter system 
consisting of scintillating-pad and preshower detectors, an electromagnetic calorimeter 
and a hadronic calorimeter. Muons are identihed by a system composed of alternating 
layers of iron and multiwire proportional chambers 1^. The online event selection is 


^Charge conjugation is implied throughout the paper. 
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performed by a trigger 26 , which consists of a hardware stage, based on information from 


the calorimeter and muon systems, followed by a software stage, which applies a full event 
reconstruction. 

In the B'^ —)■ J/-^7r+ and —)■ -^(2S')7r+ decay channels, the J/V’ and 'ip{2S) mesons 

are reconstructed through their decays into two muons. At least one muon with high px is 
required in the hardware trigger. The software trigger requires a charged particle with 
Pt > 1-7 GeV/c, or pt > 1 GeV/c if identified as a muon; alternatively a dimuon trigger 
requires two oppositely charged muons with px > 500MeV/c, and the invariant mass of 
the muon pair greater than 2.95 GeV/c^. 

Further offline selections require a good quality muon track with px > 550MeV/c, a 
good quality vertex for the reconstructed J/'0 or '0(2*5') candidate, and the reconstructed 
J/'ip and ‘ip{2S) masses to be within ±100MeV/c^ of their known values 27 . The mass 
resolution for both resonances is 14MeV/c^. The muon track pair and the pion track are 
required to be inconsistent with originating from a PV. The pion track is required to be of 
good quality, and to have a px greater than 500MeV/c. The particle identification (PID) 
information for pions is used to reduce the contamination from kaons and protons. The 
B^ candidate is required to have a good quality vertex and a reconstructed mass within 
±500MeV/c^ of its known mass 27 , which corresponds to more than ten times the mass 
resolution. To further separate signal from background a boosted decision tree (BDT) 
selection using the AdaBoost algorithm 28,29 is applied. The selection uses more input 


variables and is more sophisticated compared to the previous analysis 

Simulated samples are generated to study the behaviour of signal events. The B^ signals 
are generated with a dedicated generator Bcvegpy 30,31 through the dominant hard 
sub-process gg —)■ B^ + b + c. The fragmentation and hadronisation processes are simulated 
with P YTHIA 32,33 . The detector simulation is based on the Geant4 package 34,35 


The BDT classifier uses information on the candidate’s kinematic properties, decay length, 
vertex quality, impact parameter and angle between the particle momentum and the vector 
from the primary to the secondary vertex. The distributions of the variables that are 
used in the BDT are similar for B'^ —>■ and B'^ —>■ 'ip{2S)TT^ decays. The simulated 

sample of B'^ —)■ J/ipn^ is used as the signal sample for the BDT training. The main 
background is combinatorial, and is represented by the upper sideband in the mass 
spectrum from the 50 —>■ J/'07r+ data sample, requiring the reconstructed mass to be in 
the range [6346, 6444] MeV/c^. Since the upper sideband is used for the BDT training, the 
BDT could over-perform in this region and distort the expected combinatorial background 
in the signal region. To avoid possible bias, two BDT classifiers are trained, denoted as 
BDTl and BDT2 in the following. The 50 —)■ J/07r+ simulation and data samples are 
both split into two halves. One half of the simulated data sample and of the 50 upper 
sideband is used to train the BDTl classifier, and the other half for BDT2. Each BDT 
classifier is applied to the other half of the 50 —>■ J/xjjTT^ data sample, which is not used 
for its training. The 50 'ip{2S)Ti^ data sample is also split into two sub-samples, one 
for each BDT classifier. The threshold value for the BDT response is chosen to maximise 
the signal significance. Finally, the gT invariant mass window [3030, 3170] MeV/c^ is 
applied to J/'ifj candidates, and [3620, 3760] MeV/c^ to 0(25') candidates. 
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After the full selection, the background in the —)■ J/'07r+ sample consists of three 

categories: combinatorial background; partially reconstructed background, mainly from 
—>■ decays with —)• 7r~^7T^, where the vr® is not reconstructed; and contamina¬ 
tion from the Cabibbo-suppressed decay, B^ —)■ with the kaon misidentihed as 

a pion. The background in the B^ —)■ 'ip{2S)7T~^ sample consists of a combinatorial back¬ 
ground and a partially reconstructed background. The contribution from B^ —)■ 'ip{2S)K~^ 
is negligible. 

The signal yields are extracted from unbinned extended maximum likelihood hts to 
the invariant mass distributions of or 'ip{2S)Ti^ in the range [6027, 6527] MeV/c^, as 

shown in Figs. and for 2011 and 2012 data, and are summarized in Tables and [^ To 
improve the B'^ mass resolution, the masses of J/0 and 'ip{2S) candidates are constrained 
to their known values j^. For the B^ —)■ channel, the signal probability density 

function is modelled by the sum of two double-sided Crystal Ball functions 36 , with 


the same mean value and tail parameters determined from simulation; the combinatorial 
background is described with an exponential function; and the partially reconstructed 
background is modelled with the distribution of the B'^ invariant mass obtained from 
a simulated B^ —)■ J/^ p~^ sample using a kernel estimation 37 . This last shape is 


convolved with a Gaussian distribution to take into account a difference in mass resolution 
between data and simulation. For the B^ —)■ background, the shape of the B'^ 

mass distribution is modelled by a double-sided Crystal Ball function with parameters 
determined from simulation. For the B'^ —>■ 'ip{2S)Ti^ channel, due to the limited statistics, 
the signal shape is modelled by a single double-sided Crystal Ball function with the tail 
parameters determined from simulation; the combinatorial and partially reconstructed 
backgrounds are described with the same models as used for the B^ —)■ channel. 

The total selection efficiency is the product of the detector geometrical acceptance, 
the trigger efficiency, the reconstruction and selection efficiency, the PID efficiency, and 
the BDT classiher efficiency. All efficiencies are determined using simulated samples. To 
account for any discrepancy between data and simnlation, the PID efficiencies are calibrated 
using a tt’*' sample from Z7*-tagged —)■ K~ti^ decays. The BDT classifier efficiencies of 
BDTl and BDT2 are slightly different. After correcting for the BDT classiher efficiencies 
the signal yields of the sub-samples are consistent within the statistical uncertainties. The 
BDT classiher efficiency, ^bdt, and the product of all other efficiencies, e*, are listed in 


Table 1: Summary of the signal yields and efficiencies for the —>■ J/'07r'*“ decay. 


2011 


BDTl 


BDT2 


2012 


BDTl 


BDT2 


yield 

^BDT 


437 ± 24 475 ± 26 

(62.99 ± 0.07)% (69.29 ± 0.06)% 

(1.392 ±0.003)% 


883 ± 34 950 ± 36 

(62.33 ± 0.06)% (68.50 ± 0.06)% 

(1.339 ±0.003)% 
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Figure 1: Fit to the reconstructed mass distribution for —>• using (top) 2011 and 

(bottom) 2012 data samples. The plots on the left (right) correspond to the data selected with 
BDTl (BDT2). Black points with error bars represent the data, and the various components are 
indicated in the keys. 


Tables [T] and [U 

Several sources of systematic uncertainty on the Rq measurement are studied and are 
summarized in Table To account for the uncertainty due to the signal shape modelling, 
the data are rehtted with an alternative shape. The B'^ invariant mass distributions are 
modelled by a kernel estimation convolved with a Gaussian function, as determined from 
simulation. A difference of 0.6% from the nominal result is observed and is taken as a 
systematic uncertainty. 

The modelling of the partially reconstructed background can also introduce a systematic 
uncertainty. This is estimated by reducing the ht range to [6164, 6527] MeV/c^ to exclude 
its contribution. A change of 2.4% in the result is observed. In the nominal hts, the 
parameters for —>■ and the partially reconstructed background are hxed; the 

results change by less than 1% when these parameters are allowed to vary. The systematic 
uncertainty due to background modelling is estimated to be 2.4%. 

Systematic uncertainties on the measurement can be introduced by the BDT 
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Figure 2: Fit to the reconstructed mass distribution for 'ip{2S)7r~^ using (top) 2011 and 

(bottom) 2012 data samples. The plots on the left (right) correspond to the data selected with 
BDTl (BDT2). Black points with error bars represent the data, and the various components are 
indicated in the keys. 


classifier efficiency if the simulation fails to describe the data. The distributions of all 
training variables from simulation and background-subtracted data are compared, where 
the background subtraction is performed using the sP/ot technique, taking the invariant 
mass as the discriminating variable 38 . They are generally in agreement within statistical 


fluctuations. Only one variable, which describes the consistency between the pion track 
and the PV, indicates small differences between simulation and data. Therefore, the 
simulated sample is reweighed to match the data, and the BDT efficiencies are recalculated 
with the reweighed simulated sample. The result obtained with these BDT efficiencies is 
different from the nominal value by 0.2%, which is taken as the uncertainty from the BDT 
classifier. 

The efficiences determined from simulated samples have uncertainties due to the limited 
statistics. This leads to an uncertainty of 0.3%. An uncertainty of 1.1% is assigned due to 


imperfect simulation of the trigger, which is determined using data driven methods 39,40 
The B'^ lifetime of simulated samples is set according to the latest LHCb measurement^ 
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Table 2: Summary of the signal yields and efficiencies for the B'^ —)• decay. 


2011 


2012 


yield 

^BDT 


BDTl 


BDT2 


14.4 ±4.5 19.6 ±5.3 

(58.79 ±0.11)% (65.84 ±0.11)% 
(1.631 ±0.006)% 


BDTl 


BDT2 


40.1 ±7.1 30.8 ±7.0 

(58.32 ± 0.08)% (65.08 ± 0.08)% 

(1.529 ±0.005)% 


To estimate the systematic uncertainty due to this, the B'^ lifetime is varied within 
the uncertainty of this measurement, and the change in the result, 0.1%, is taken as a 
systematic uncertainty. The total systematic nncertainty is 2.7%. 

The ratio of the branching fractions with J/'0 and '0(25) mesons decaying to dimuons, 
denoted as 

B{Bt ^ i,{2S),T*,i,(2S)^ 

B{B+ —)■ J/07r+, J/0 —)■ /r+/i“) 

is calcnlated as 

„ _ iV|S;i(g+ ^ ii2S)n*) + iva,(g+ ^ ii2S)n*) 

-> J/^ir+) + J/i/>x+) 

where ^ 2011 ( 2012 ) signal yields from 2011(2012) after efficiency correction. The 

ratio is measnred to be 


R = 0.0354 ± 0.0042 (stat) ± 0.0010 (syst). 

The ratio of the branching fractions of —)■ '0(25)7r’'' and —)■ is calculated as 

BiJ/i) /r+/i“) 


Rk — R X 


B{'ip{2S) —)■ /r+/r' 


( 4 ) 


Assuming electroweak universality, the J/0 —?• and "0(25) —)■ branching 

fractions can be substituted with the more precisely measured ones in the e’''e“ channel 27 . 
Using these valnes, the ratio Rb is measnred to be: 

Rb = 0.268 ± 0.032 (stat) ± 0.007 (syst) ± 0.006 (BF), 

where the hrst uncertainty is statistical, the second is systematic, and the last term is 
due to the uncertainty on -B(J/'0 —)■ e'^e~)j B(^{2S) —)■ e+e“). This result is in agreement 
with the previous LHCb result j^. Our measurement is consistent with the predictions 
of non-relativistic QCD at next-to-leading order (0.26lo;og) 18 and perturbative QCD 
based on factorization (0.29 ](q;}[) |1^. The resnlt disfavors the theoretical calculations 
based on the relativistic qnark model 1141, the qnark potential model 15 , the relativistic 


constitnent qnark model 16 and the QCD relativistic potential model 17 


6 





















Table 3: Summary of systematic uncertainties on the Rq measurement. 


Component 

Uncertainty 

Signal shape 

0.6% 

Background shape 

2.4% 

BDT classiher 

0.2% 

Monte-Carlo Statistics 

0.3% 

Trigger efficiency 

1.1% 

Bf lifetime 

0.1% 

Total 

2.7% 
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